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The Goldstone theorem states that there should be a massless mode for each spontaneously 
broken symmetry generator. There is no such rotational mode in crystals, however superconducting 
quantum nematics should carry rotational Goldstone modes. By generalization of thermal 2D defect 
mediated melting theory into a 2+1D quantum duality, the emergence of the rotational mode at 
the quantum phase transition from the solid to the nematic arises as a deconfinement phenomenon, 
with the unusual property that the stiffness of the rotational mode originates entirely in the dual 
dislocation condensate. 



PACS numbers: 11.30.Qc,61.30.Gd,m74.25.N- 

Introduction. Crystals break both translations and 
rotations but its Nambu-Goldstone modes (phonons) are 
translational only. Where are the modes expected from 
the Goldstone theorem P~H3] associated with the rota- 
tional symmetry breaking [3HS]? Their absence is crucial 
to engineering since the presence of massless 'rotational 
phonons' would imply an elastic response towards torque 
stresses, corrupting crankshafts and so forth. However, 
matter can also just break space rotations leaving trans- 
lations intact. Such nematic liquid crystals [5J [7J [8] do 
exhibit an elastic response towards torque stress. In reg- 
ular nematics this issue is muddled by the peculiarity 
that the rotational Goldstone modes are overdamped due 
to the decay in hydrodynamical rotational flows [SI [H]- 
In a recent development evidences for the existence of 
zero-temperature quantum liquid crystals have been ac- 
cumulating [10l [11] . A particular variety appears to be 
present in cuprates [I2HI5] and iron-arsenides [HI [TjJ 
with ground states that also superconduct. Given that 
the hydrodynamical vorticity is gapped in superconduc- 
tors, this is a natural theater to look for propagating 
rotational Goldstones. 

The conventional theory of liquid crystals departs from 
a kinetic gas theory perspective [7j [8], a view which is 
also popular in the quantum realms [111 I18j . However, 
a dense nematic might be closer to a system which is 
locally still like a solid, turning at long distances into a 
fluid because dislocations (translational topological exci- 
tations) have proliferated. This notion was introduced by 
the famous Kosterlitz-Thouless-Nelson-Halperin- Young 
theory of the thermal melting of a 2D floating crystal 
into the hexatic liquid crystal |19 21j. More recently, 
a powerful field theoretical weak-strong duality [2"2Tf2~l] 
was mobilized to address the zero-temperature quantum 
melting of a bosonic solid into a nematic superfluid [25l - 
|2"8] . This follows closely the pattern of the superfluid 



superconductor Abelian-Higgs duality in 2+1D [25l I29j . 
The rigidity of the crystal phase is captured by gauge 
fields (here called stress photons) representing the ca- 
pacity of the medium to propagate shear forces. These 
are sourced by the dislocations, and the quantum melt- 
ing into the nematic superconductor corresponds to a 
proliferation and Bose-condensation of dislocations in an 
effectively relativistic 'stress superconductor' explained 
below, which turns out to be an electromagnetic super- 
conductor as well [23 1 [27 1 [25], 

This theory can be regarded as universal for the long- 
wavelength regime in the adiabatic continuity sense. 
Here we want to present the counterintuitive but elegant 
resolution it offers for the phonon conundrum of the first 
paragraph. We will demonstrate that the 'crankshaft 
rigidity' of the solid is explained as a confinement of 
torque stress, in literally the same guise as color forces 
are confined in the IR of a non-Abelian gauge theory. At 
the quantum phase transition to the superfluid nematic, 
the shear forces acquire a Higgs mass 51 signaling the dis- 
appearance of the phonons due to the restoration of the 
translations. The associated length is the 'shear penetra- 
tion depth' having the same status as the London pene- 
tration depth but now specifying the length over which 
shear forces decay in the quantum nematic, in turn equiv- 
alent to the correlation length. But at the same time, 
the dislocation condensate deconfines the torque stress 
indicating the appearance of a massless rotational Gold- 
stone mode, going hand-in-hand with the deconfinement 
of the sources of torque stress: the disclinations (rota- 
tional topological defects). 

A quite non-trivial prediction follows from the duality: 
the quantities specifying the collective excitations of the 
superconductor coincide with those describing the crys- 
tal elasticity [Eq. ([IJ] and in addition only the shear 
Higgs mass is needed. The velocity of the rotational 
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Goldstone is set by the phonon velocity but the rota- 
tional (nematic) modulus diverges with Q~ 2 approach- 
ing the quantum phase transition back to the solid, see 
Eq. (jTsj) . The mechanism explaining this single scale 
is remarkable: it is rooted in a new form of 'Higgsing'. 
The rotational stiffness in the quantum nematic is actu- 
ally carried by the dislocation condensate — as opposed 
to as a remnant of the solid — knowing about the ro- 
tational symmetry breaking through the orientations of 
the dislocations. The 'longitudinal Goldstones' (gauged 
phase modes) of this special condensate rearrange in such 
a way that, besides 'being eaten by the stress photon', 
they also form the massless rotational Goldstone mode. 
Let us use the remainder to explain how this works. 

Dual elasticity. We depart from the minimal theory 
of quantum elasticity, describing an isotropic medium in 
2+ ID. One combines the gradients of the displacement 
field u a (x) in symmetric strain tensors w a b = \{d a u b + 
dbu a ). The Lagrangian describing the static elasticity 
as well as the acoustic phonons includes a kinetic term 
(t = it is imaginary time) |25j . 

C = \ P {d T u a ) 2 + ^KJ 2 + ~(w ab w ba ), (1) 

in terms of the compression modulus k, the shear mod- 
ulus p and the mass density p, such that the longitu- 
dinal and transverse phonons propagate with velocities 
c£ = (ft + p and c\ = p/p, respectively. To establish 
contact with gauge field theory one focuses on the capac- 
ity of the medium to propagate forces and this is accom- 
plished by strain-stress duality. By Legendre transfor- 
mation one defines the stress tensor af, = x ,i c ^ , where 

= (-^d T ,d m ). The dual stress field satisfies the con- 
servation law and Ehrcnfest constraint, 

= °- ( 2 ) 

tarn&m = 0- (3) 

The Ehrenfest constraint Eq. ^ eliminates spontaneous 
torque, in dual correspondence with the absence of anti- 
symmetric 'rotational strains' u> — ^(d x u v — d y u x ) in 
the leading gradient terms defining the action Eq. ([!]). 
Although less familiar in this context, one can introduce 
gauge fields in the same way as in Abelian-Higgs duality 
by exploiting the fact that the conservation of stress Eq. 
([2]) can be imposed in 2+1D by parameterizing the stress 
tensor in terms of 'flavored' (spatial labels a) U(l) gauge 
fields b a x as [22], 



<(x)=e fl vxd v bl(x). 



(4) 



In the case of Abelian-Higgs (or XY) duality the unique 
internal sources for the dual photons are the vortices, 
with the implication that the 2+ ID superfluid is mapped 
onto the electromagnetic Coulomb phase. Similarly, here 
one finds that the stress photons are sourced by disloca- 
tion currents: £ sourcc = 6° J£, where J° = e^xd^dxu^ 



and Wg ing is the singular part of the displacement field 
[22l [24]. Compared to the vortices, this structure is richer 
since the dislocations carry Burgers vectors as topologi- 
cal invariants. These are enumerated through the flavor 
labels a, encoding that dislocations only restore transla- 
tions in the direction of their Burgers vectors. The bot- 
tom line is that phonons can actually be viewed as the 
duals of stress photons, describing the exchange of elastic 
forces between the translational topological defects. 

Torque stress gauge field. The breaking of rotational 
symmetry is completely implicit in the above derivation. 
Its influence is embodied by the disclination, the topo- 
logical defect uniquely associated with the restoration of 
the isotropy of space. Our main discovery is a general- 
ization of stress gauge field theory, rendering the role of 
rotational symmetry explicit in terms of a torque stress 
gauge field. For this we have to treat the Ehrcnfest con- 
dition Eq. ([3| as a dynamical constraint. This imposes 
that only the symmetric part cr^ = o™ of the full stress 
tensor = af n + <7^ is present, while its antisymmetric 
part (jj^ — ~®""a vanishes. Our trick is to use the an- 
tisymmetric strain ui as an independent dynamical field 
by implementing the relation uj — -^CijdjU^ as a con- 
straint. We conveniently choose CmaO'm as the Lagrange 
multiplier and add the term e ma (jj^(u; — ^CijdiU^) to the 
Lagrangian Eq. dualizing into 



C 



constraint 



(5) 



As usual the conservation of stress Eq.Q is imposed by 
u a as a Lagrange multiplier [5SJ [38] , but now also the 
Ehrenfest constraint Eq. (|3| is incorporated by integrat- 
ing out the smooth part of w. Now comes the crucial step: 
substitute Eq. Q into Eq. Q to obtain the conserved 
torque stress r^, 



d^T,,. = 0. 



(6) 
(7) 



The conservation of torque Eq.Q can be imposed by 
parameterizing it in terms of a vector U(l) torque stress 
gauge field associated with the the 0(2) rotational 
symmetry of 2D space, 



Tu = £uv\d v h\. 



(8) 



Substituting the dual gauge fields Eq. Q and Eq. Q 
into Eq. |5]), we obtain for the full dual Lagrangian, 

r ^ ^o, a 1 / a a i a a b\ 

L = ~Yp r r ~ m m + m ° - 

+ b a x J a x + h x Q x + < m (d^) + w sm (e am <7° ) (9) 

This is just the standard action for elasticity written 
in stress fields, having the constraints Eq. ([2]), ^ 
explicit in terms of the Lagrange multipliers Ug m ,u; sm . 
These are resolved by the stress gauge fields where b x 
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is sourced as usual by the dislocation current, but the 
novelty is the term representing the disclination current 
@a = £\ t j,i>d l j,dvUj s i ng , enumerating the singular part of 
the rotation field [22, 24 sourcing the torque stress gauge 
field h\. The physical stress fields <r^ and r M are invari- 
ant under the following gauge transformations involving 
arbitrary fields f a and g, 



b a x + d x f a 



(10) 
(11) 



Substituting into Eq. |9L we find the conservation laws 
for the defect currents [22] , 



0. 



(12) 
(13) 



which appear here as natural consequences of the gauge 



symmetries Eq. (10l-(ll). 

It is an old wisdom that disclinations are infinite- 
energy defects in solids (the crankshaft effect), but we 
are now in the position to address this in gauge the- 
ory language. After substituting Eq. Q,([9| acquires 
a Maxwell-like form. Going to momentum space where 
p — (— w,q), by exploiting Eq. ([6]) in combination with 
Eq. ^ one straightforwardly finds that the part of the 
Lagrangian involving torque stress gauge fields has the 
form c^h\h T + h T Q T . This implies that for instance 
the static interactions between disclinations behave as, 



(h T ( P )h T (0)) 



1 



5 2 Z[Q] 



Z 60 T {Q)6@ T (p) 



■H 1 



2k/j, q 4 



(14) 



The energy of a static disclination-antidisclination pair 
therefore increases quadratically with their separation. 
This is quite like the confinement phenomenon in non- 
Abelian Yang-Mills theories, except that the potential 
there is linear. Although the mechanism is in detail 
rather different, the physical meaning of the confinement 
phenomenon is very similar in both cases. In QCD the 
quarks are the sources of color forces mediated by glu- 
ons, but the confining state is 'color-undeformable', only 
allowing for color singlets, while quarks and gluons do 
not have physical existence in the infrared. In the solid, 
disclinations are like quarks and the torque gauge fields 
are like gluons but the vacuum of the solid is a 'torque- 
gauge singlet' and can therefore be used as crankshaft. 

Quantum nematic. Can the torque gauge field be de- 
confined? The reader might anticipate the answer: in the 
superconducting quantum nematic, disclinations should 
occur as the topological defects, while our torque stress 
gauge field should dualize into the massless rotational 
Goldstonc mode. We will show this explicitly mobilizing 
the Abelian-Higgs duality technology which directly ap- 
plies to the problem of quantum melting of the bosonic 
crystal into the nematic superfluid in 2+ ID. This is in 



essence just the quantum version of the 2D topological 
melting theory developed in the 1970s 19 21 . Dislo- 
cations are uniquely associated with the restoration of 
the translation symmetry. The nematic fluid can there- 
fore always be considered as a crystal where the dislo- 
cations have spontaneously proliferated. In 2+1D the 
dislocations are themselves bosonic point particles and 
when they proliferate they will Bose-condense. Since the 
stress gauge fields mediating long-range interactions be- 
tween the dislocations are like photons, this dislocation 
condensate is a (relativistic) 'dual stress superconduc- 
tor'. Elsewhere we have explored the gross properties of 
such bosonic quantum liquid crystals 25, 29, 30 . In short 
summary, the shear photons acquire a Higgs mass f2, indi- 
cating that the shear forces in the nematic become short 
ranged after restoring translational symmetry, while in 
principle a massive, propagating shear photon should be 
present [26] [27j [30] . Since compressional stress (as op- 
posed to shear stress) does not couple to dislocations [31] , 
the longitudinal phonon of the solid turns into a pure 
zero-sound mode, that coincides with the phase mode 
of the superfluid. Upon coupling to electromagnetism it 
can also be shown that the dislocation condensate shows 
a genuine Meissner effect — it is a regular superconduc- 
tor [251 EH [30]. 

Lacking the torque gauge field formalism introduced 
in the above, in this older work the fate of the rota- 
tions/torque in the dislocation condensate was rather ob- 
scure. Using torque gauge fields this now becomes explic- 
itly tractable and we find that the mechanism by which 
the rotational Goldstone is 'born' in the nematic state is 
elegant but counterintuitive: it is formed from the longi- 
tudinal modes of the dislocation condensate itself. This is 
rooted in two special properties of this condensate. First, 
as in all Abelian-Higgs type dualities, there is only one 
velocity (c.q. the phonon velocity) and such condensates 
are truly relativistic Higgs phases where the condensate 
phase mode propagates with the same velocity as the 
photons. Therefore, the 'condensate mode is eaten by 
the gauge field', giving rise to an additional (massive) 
'longitudinal photon'. Secondly, the dislocations carry 
Burgers vectors and therefore the condensate has multi- 
ple components. The rotational point group symmetry 
imposes a unique way to construct the dislocation con- 
densate corresponding with the nematic [26-28, 32l I33j : 
the dislocation Burgers vectors have to orient along the 
directions associated with the crystal point group and 
the nematic condensate is formed by populating all these 
directions equally. Taking these requirements it follows 
from symmetry reasoning that the Josephson form of the 
Higgs term of the nematic condensate has to have the 
form, 
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The dislocation condensate fields are $ a 
where a indicates the Burgers directions, £1 - 



(15) 

= |$o| e i&> 
|$ a | is the 
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Higgs mass and <j) a is the phase variable of the condensate 
field. The complete action for the quantum nematic is 
obtained by adding Eq. Q and omitting the dislocation 
source term. Focusing on the symmetric gauge and con- 
densate fields (^symmetric) one discovers the physics dis- 
cussed in the previous section. However, with the torque 
stress fields at hand it becomes easy to address the rota- 
tional sector as well. The novelty is that the Higgs term 
Eq. (151, which usually encodes for mass, is now the 



origin of the massless rotational Goldstone mode! 

The proof is easy. By taking the unitary gauge fix 
b\ — > + d\4> a the condensate modes are shuffled into 
the gauge sector and the Higgs term becomes, 



C 



Higgs 



2/i 



(16) 



According to Eq. ^ the torque stress can be expressed 
in terms of the 6" as t t = — b% — b^.,T x = b% and 
T y = b^.. The full action can be rewritten as £f u n = 
^rotational + ^symmetric ■ The long- wavelength rotational 
dynamics is just governed by the Higgs term Eq. ( 16 ) 
acquiring eventually the very simple form, 



C 



rotational 



1 n 2 ~ - 



(17) 



where F,. 



d u h^ is parameterizing the torque 
stress term t^t^ in terms of the torque stress gauge fields 
being sourced by the disclination currents 0^. As a 
caveat, an extra factor 1/2 showing up in front of the t t t t 
terms causes the velocity of the rotational Goldstone to 
be reduced relative to the transv ersa l phonon velocity by 
cr = Ct/v2- The form of Eq. (17), besides periodicity 
issues that are kept implicit here, indicates that disclina- 
tions in 2+1D are just like vortices, and the rotational 
Goldstone mode has the same role as the phase mode 
in the superfluid, giving rise to the usual logarithmic in- 
teractions. However, the origin of the rotational rigidity 
is remarkable. Dualizing backward to the 'strain side', 
the elastic action associated with the rotational rigidity 
becomes that of a free massless scalar field 6, 
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-d T (j>) 2 + (d m cj>f 



(18) 



the rotational mode is encoded in the phase modes of the 
dislocation condensate itself. This in turn makes sense 
since the dislocation condensate communicates via its 
Burgers vectors with the rotational symmetry breaking 
and it should therefore sustain a restoring force against 
a rotational twist. 

Outlook. In the seemingly mundane context of elas- 
ticity of crystals and liquid crystals we have uncovered a 
remarkably rich and elegant gauge field-theoretical struc- 
ture, explaining the rotational rigidity of the superfluid 
liquid crystal in terms of a deconfinement transition from 
the confining solid. How to test our quantitative predic- 
tion Eq. ( 18 ) experimentally? The available supercon- 



ducting nematics [THUS, -^J may be of the strongly cor- 
related kind, but the presence of an ionic lattice breaking 
the space rotations explicitly is an obscuring factor. This 
'lattice pinning' of the (uniaxial) director will cause an 
anisotropy gap in the spectrum of the rotational mode, 
which has to be small compared to the other scales in or- 
der for the continuum theory to be of relevance. A more 
fundamental complication is that one can only probe 
the electron nematics through their electromagnetic re- 
sponse. The massive shear mode is observable in linear 
response [27] , albeit in a kinematical regime which is be- 
yond the reach of present spectroscopic techniques. Sim- 
ilarly, it is straightforward to compute [551 HZ] the cou- 
pling of external electromagnetic fields to the rotational 
mode: it is in principle visible in the transverse dielec- 
tric function with the complication that h\ only cou- 
ples to the gradient of the magnetic field B as h^d^B, 
where T and L are the transverse and longitudinal di- 
rections. This gradient coupling greatly complicates the 
direct measurement of this dynamical mode and we leave 
it as a challenge for the community to find out how to 
create and interrogate the 'platonic' nematic supcrfluids 
of the weak-strong duality in the laboratory. 
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The conclusion is that besides the shear modulus /i and 
the mass density p of the background crystal (recall that 
c r = c|/2 = p/2p) this just depends on the Higgs mass 
O of the dislocation condensate. The rotational modulus 
K = 2fi/Q 2 is diverging upon approaching the quantum 
phase transition to the crystal, which is perfectly reason- 
able given that the crystal is completely undeformable in 
this regard. 

Although the derivation is simple, the physics it de- 
scribes is quite complicated. The mathematics shows 
that the rotational rigidity of the quantum nematic is 
governed by the Higgs term. As one can explicitly check, 
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